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ABSTRACT 

Recently the Spitzer Space Telescope observed the transiting extrasolar planets, TrES-1 and 
HD209458b. These observations have provided the first estimates of the day side thermal flux from two 
extrasolar planets orbiting Sun-like stars. In this paper, synthetic spectra from atmospheric models 
are compared to these observations. The day-night temperature difference is explored and phase- 
dependent flux densities are predicted for both planets. For HD209458b and TrES-1, models with 
significant day-to-night energy redistribution are required to reproduce the observations. However, 
the observational error bars are large and a range of models remains viable. 
Subject headings: Planets: exoplanets, radiative transfer 



1. INTRODUCTION 

Of the more than 100 extrasolar planetary systems dis- 
covered so far, only 7 have near edge-on orbits. These 
transiting planets are crucial for understanding giant 
planets in general since their masses and radii can eas- 
ily be determined, and careful multi- wavelength observa- 
tions can r eveal some information about the planet's at- 
mosphere llBrown et all 120011: iCharbonneau et all 120021: 
iVidal-Madiar et alJ200l ~ Recently, two transiting plan- 
ets, HD209458b and TrES-1, were observed with the 
Spitzer Space Telescope, providing the first direct mea- 
surements of their thermal flux. By comparing IR fluxes 
in and out of secondary eclipse (when th e planet is be- 
hind the star 'I . ICharbonneau et al.l l)2005l) measured the 
planet-star flux density ratio a t 4.5 a nd 8 /jm for TrES-1. 
Independently, iDeming et alJ i|2005T) measured the flux 
density ratio at 24 fim for HD209458b. These measure- 
ments provide the best constraints, so far, on the thermal 
structure and chemical composition of highly irradiated 
EGPs. 

HD209458b and TrES-1 both have very short peri- 
ods (just a few days) and orbital separa tions less than 
0.05AU ijHenrv et alJl2000l lAkmso et alJl2004j) . At such 
small orbital separations, they are substantially heated 
by radiation from their parent stars. Following the dis- 
cover y of the planet around 51 Pegasi ijMavor fe Quelozl 
1995), a variety of atmosphere models suitable for esti- 
mating the pro perties of these so-called "hot Jupiters" 
were published ISeaeer fc Sassek~>vlll998t iBarman et al.l 
I200T1 ISndarskv et alfl l2003L However, since these plan- 
ets probably have strong day-to-night photosphcric dif- 
ferences, their potential lack of symmetry adds compli- 
cations to an already difficult model atmosphere prob- 
lem. Several of the most challenging issues are the cov- 



erage and types of clouds (if present), redistribution of 
the absorbed stellar flux by atmospheric currents, depth 
dependent non-solar abundances, and photospheric tem- 
perature and pressure gradients from the day to night 
sides. These problems have been dealt with (or avoided) 
in a variety of ways and, therefore, a variety of model 
predictions exist. 

New models, that estimate the horizontal atmospheric 
gradients under the assumption of radiative-convective 
equilibrium, are presented below. Several of the stan- 
dard assumptions for global energy redistribution are ex- 
plored. Model results are also compared to the most re- 
cent Spitzer measurements and estimates are given for 
the planet-star flux density ratios in the Spitzer 24 /urn 
MIPS band and the four IRAC bands as a function of 
orbital phase. 

2. MODEL CONSTRUCTION 

The irradiated models presented below were 
calculated using the PHOENIX atmosphe re code 
(|Hauschildt fc Baronl fToOflt lAllard et al] 12001ft adapted 
to include extrinsic radiation as described in Barman et 
al. (2001; here after BHA01, 2002) . The spherically 
symmetric radiative transfer and chemical equilibrium 
equations were solved self-consistently, while explicitly 
accounting for the wavelength-dependent extrinsic 
radiation. The extrinsic radiation was also modeled 
with PHOENIX and, in each case, closely matches the 
observed parent star spectrum. 

The major differences between the computation of the 
BHA01 models and those presented here lie with the 
treatment of dust in the atmosphere and the assump- 
tions concerning the redistribution of absorbed stellar 
flux over the planet's day and night hemispheres. 
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Table 1. Model Parameters 



HD209458bt TrES-1* 



T* (K) 
R* (R ) 

M* (Mq) 


6088 ± 56 
1.145 ±0.049 
1.06 


5250 ± 75 
0.83 ± 0.05 
0.89 


Rp (Rjup) 

M p (M Jup ) ■ 
a (AU) 


i 49+0- 10 

0.69 
0.0468 


0.76 
0.0393 



Ref ere nces. ( f) IRibas et alJ 120031) . ICodv k. Sasselovl 

dOOD; {%) ISozzetti et all 120041) 



2.1. Cloud-free Assumption 

Unlike many brown dwarfs, EGPs do not necessar- 
ily have convective photospheres. For EGPs with small 
orbital separations, irradiation can suppress convection 
to depths well below the photosphere, leading to a 
fully radiative atmo sphere across most of the day side 
lIGuillot et al.l l"l996). Since radiative photospheres have 
short sedimentation time-scales, cloud growth should be 
difficult to sustain making the cloud-free assumption rea- 
sonable. However, this assumption may break down 
if strong zonal winds are present and advective time- 
scales are comparable to sedimentation time-scales. The 
efficiency of gravitational settling will also de pend on 
the p oorly constrained eddy diffusion coefficient ijRossowl 
1978). 

In BHA01, cloud- free atmos pheres were m odeled us- 
ing the "Cond" opacity setup l|Allard et all 120011) . The 
Cond setup accounts for dust formation in the atmo- 
sphere, as determined by chemical equilibrium equations, 
but excludes the dust opacity when computing the ther- 
mal and spectroscopic properties of the atmosphere. Ex- 
cluding the dust opacity was intended to approximate 
the effects of cloud formation followed by efficient grav- 
itational settling (often referred to as "Rainout") which 
acts to deplete an atmosphere of many important re- 
fractory elements. Rainout and the depletion of refrac- 
tory elements have been recognized for some time as 
important processes in the atmospheres of Jovian plan- 
ets and brown dwa rfs ( Fcel ev fc LodderslllQQfil: iLoddersl 
119991 Burrows et al.ll2000UMarlev et alJl2002|) . The suc- 
cess of these models is strong motivation for exploring the 
rainout assumption in extrasolar planet atmospheres. 

While the Cond approximation does effectively remove 
refractory elements, it does not alter the overall abun- 
dance of an element sequestered by grain formation at 
a given temperature and pressure. Grain formation and 
efficient gravitational settling alter the abundance of an 
element by continually removing the refractory elements 
that make up a grain until the grain is no longer able to 
form due to a lack of one or more constituents. In the 
absence of replenishment (e.g., by convective updrafts), 
certain layers of an atmosphere where only gravitational 
settling has occurred will not only be free of dust species 
but also free of a significant number of metals and related 
molecules. Fo r a description o f this effect in Jupiter's at- 
mosphere, see ILoddersl i|1999l) . 

The Cond case, therefore, does not go far enough in re- 
moving refractory elements from the upper atmosphere 




Fig. 1. — Illustration of the planet's hemisphere divided into 
small regions, each to be modeled independently. Each concentric 
band around the sub-stellar point (depicted here as the convergence 
point of longitudinal lines) receives the same amount of incident 
stellar flux at the same incident angle. Therefore, each band is 
assumed to have the same atmospheric structure and emergent 
intensities. The angle between the incident ray and the outward 
pointing surface normal is 9. When computing the emergent flux, 
only those intensities aligned with the observer line-of-sight are 
integrated over the visible hemisphere. The angle between the 
surface normal and an observer's line-of-sight is (3. 



and can lead to a small concentration of key absorbers 
like TiO and VO that are important in irradiated at- 
mospheres. These two molecules have strong absorption 
bands near the peak flux densities of solar type stars. 
Consequently, their presence or absence can greatly af- 
fect the depth at which the stellar flux is absorbed in 
the planet's atmosphe re and, thus, alter t he predicted 
atmospheric structure ijHubenv et al.ll2003|) . 

In this work, an improved cloud-free model was used 
that iteratively reduces (at a given layer) the elemen- 
tal abundances involved in grain formation and recom- 
putes the chemical equilibrium with each new set of strat- 
ified elemental abund ances. This mode l is similar to 
the Rainout model of iBurrows fc Sharpl l)1999|) . except 
that the depletion of elements is continued until grains 
(and thus grain op acities) are no longer present (see also 
lAllard et alJ (2003) for more details). The resulting equi- 
librium chemistry and opacity sampling of this cloud-frcc 
model are fully self-consistent, unlike the earlier Cond 
models which simply excluded the grain opacity. Also, 
for the models described below, Ti and V were signif- 
icantly depleted from the photosphere by the rainout 
process leading to negligible concentrations of TiO and 
VO. For a detailed discussion of the differences between 
Rainout and Cond irradiated models, see Barman et al. 
(2005, in preparation). 

2.2. Modeling the Day-Night Gradients 

The majority of static atmosphere models arrive at 
a single temperature-pressure (T-P) profile intended to 
represent a n average of either the day side or over the en- 
tire planet llSeager fc Sasselovll998l:IGoukenleuque et alJ 
l2TJ0fl iBaTman et al.l 120011 iSudarskv et alJl2003|) . How- 
ever, if most short period planets are well represented 
by a static atmosphere in radiative-convective equilib- 
rium, then one should expect them to have a horizontal 
(day- to- night) temperature gradient - simply due to the 
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center-to-limb variation in the amount of incident stellar 
flux received by the planet. In the absence of a 3-D model 
atmosphere code, one approach that comes closer to the 
real solution is to divide a planet's day side into a series 
of concentric regions around the sub-stellar point. In the 
static case, suitably large regions should interact very lit- 
tle via radiative transfer processes, except perhaps near 
the terminator. 

The planetary atmospheres described below were mod- 
eled by dividing the day side into 10 concentric regions 
defined by /i = cos(0), where 9 is the angle between the 
surface normal and the direction to the star (see Fig. 
For these regions, fi ranged from 1.0 (at the sub-stellar 
point) to 0.1 (the model region closest to the termina- 
tor), in steps of A/i = 0.1. The corresponding T-P pro- 
files and emergent intensities were modeled using 1-D, 
spherically symmetric, atmospheres each receiving inci- 
dent stellar flux along the appropriate angle for a given 
region. For these day side models, the radiative transfer 
equation was modified so that the incident specific in- 
tensities along any /z and azimuthal angle <fi were given 

by 

Iinc,\{^, 0) = Io,\6(jl - n')S((f> ~ <(>'), (1) 

with 5 being the Dirac delta function. In which case, it 
follows that the incident fluxes are simply, 

F inCt x(fJ>) = ^h,x = M F *' x ' ( 2 ) 

where F+^\ are the monochromatic fluxes from the star's 
surface, is the stellar radius, and d is the distance 
from the stellar surface to the planet's atmosphere. For 
the night side, a single, non-irradiated, model was used. 

All models were solved self-consistently so that each [i- 
region had a chemistry characterized by its T-P profile. 
By having chemical equilibria consistent with the T-P 
profiles across the planet's atmosphere, this approach 
naturally leads to variations in the important photo- 
spheric opacity sources from the day to night side - an 
important aspect when computing the synthetic spectra. 

Since EGPs are believed to have fully convective in- 
teriors, the intrinsic effective temperature (Tj nt ) 1 for 
each modeled region was adjusted so that, after conver- 
gence, all T-P profiles reached the same adiabat below 
the photosphere 2 . The adiabat was selected based on 
planetary interior and evolution calculatio ns for a given 
mass , age, metallicity, and irr adiation ijBaraffe et alJ 
120031 12004 iChabrier et afll2004j) . This entropy match- 
ing technique has also been used for irradiated binary 
stars and allows one to assign models with different in- 
trinsic lu minosities to different regions of the same star 
or olanet llVaz fc NordlundllT98l iNordhmd fc Vazlll99(l 
Bar man et al.ll2004j) . 

The monochromatic fluxes from the model planet were 
obtained by integrating the emergent intensities along an 
observer's line-of-sight (los) for a given observer-planet- 

In the present work, Ti nt characterizes the intrinsic luminosity 
of an irradiated model atmosphere, defined by 4irRpoTf nt . For 
non-irradiated models, the normal T e g is used to described the 
flux and luminosity. 

2 The region referred to as the photosphere lies roughly between 
P = 0.01 and 1 bar, corresponding to where the optical depth at 
IR wavelengths is near unity. 



star orientation. 

Fi os ,x = f h(9,<t>)dn. (3) 

los 

The integration was performed by distributing ~ 2000 
points over the entire surface with corresponding cuba- 
ture integration weights (wi). The distribution of points 
was determined by the minimization of potential en- 
ergy on the unit sphere. This distribution is nearly 
orientation independent, unlike the standard latitude- 
longit ude grid, and results in very small integration er - 
rors iStcinack er et alJl!996t ISloan fc Womerslevl 12001'). 
Each point on the observer-facing hemisphere was as- 
signed an emergent intensity spectrum corresponding to 
the angle, (3, between the observer's los and the surface 
normal (see Fig. [Q, and depending on the day or night 
region in which the point belonged. Numerical integra- 
tion of Eq. |3 becomes a sum over the visible points, 

Fi os ,\ = ^2 WitHhifJ'i)- (4) 

i(visible) 

In Eq. U /i = cos(0) and is not to be confused with the 
10 fi values used to divide the day hemisphere. With this 
technique, the phase-dependent spectra can be estimated 
while taking into account the center-to-limb variation of 
the planet's emergent intensities - which might include 
a combination of limb brightening and darkening. The 
emergent intensities were sampled over 114 directions per 
/x-region and were assumed to be azimuthally symmetric 
about the surface normal. Assuming azimuthal symme- 
try is justified since the present work is concerned pri- 
marily with the thermal flux, not scattered light. Phase- 
dependent optical spectra will be explored in a later pa- 
per. 

The simple approach outlined above has several limi- 
tations. Since radiation passing through one region into 
another is neglected, heating of the upper atmosphere 
near the terminator may be underestimated at low gas 
pressures. However, this is unlikely to affect the emer- 
gent thermal IR spectrum that forms deeper in the at- 
mosphere. It has also been assumed that the planet is 
tidally locked and, hence, presents a constant face to 
the parent star. The models are also time-independent 
and static and, consequently, neglect the effects of zonal 
winds that could change the thermal profiles by coupling 
the hot day side to the cooler night side. 

2.3. Energy Redistribution 

The approach outlined in the previous section is de- 
signed to model the planet's atmosphere under the as- 
sumption that the gas is truly static and in radiative- 
convective equilibrium. As such, this approach will pre- 
dict the maximum heating at the sub-stellar point and 
very little heating at the terminator. However, as men- 
tioned above, an important consequence of stellar heat- 
ing is horizontal atmospheric flows capable of transport- 
ing appreciable amounts of energy to the night side. 
The impacts of horizontal motion on EGP atmospheres 
have been modeled by a variety of groups, each predict- 
ing some level of atmospheric circulation that depends 
strongly on a dopted opacities and general approaches to 
the problem llShowman fc Guillotll2002t iCho et alJl2003t 
iBurkert et alJl2005t iCooper fc Showmanll2005j) "Despite 
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Fig. 2. — Temperature versus pressure for a sequence of irradi- 
ated atmospheres. For each model, the parameters for TrES-1 from 
Tablefflwere chosen and only the direction of the incident flux rel- 
ative to the surface normal was varied. From top to bottom, the 
models have fi = 1.0 to fi = 0.1 in steps of 0.1. The dashed lines 
indicate the approximate condensation curves for three common 
grain species. The dotted line indicates where gaseous CO and CH4 
concentrations are equal (CO is dominant to the left of this line). 
The approximate regions represented by the collection of T-P pro- 
files are shown as solid black lines on the illustrative sphere. The 
top-most T-P profile corresponds to the sub-stellar point (black 
dot on the sphere). The terminator and night side (black hemi- 
sphere) are modeled with the non-irradiated profile (lowest T-P 
curve). The radiative-convective boundary at the sub-stellar point 
and on the night side are labeled with filled circles. The thick, 
grey, dashed lines are T-P profiles for a = 0.5 (top) and a = 0.25 
(bottom) models. 



the differences in methods and results, the general con- 
sensus from these hydrodynamic simulations is that cir- 
culations can redistribute a fraction of the incident en- 
ergy over large portions of a strongly irradiated planetary 
atmosphere. 

For single, 1-D, model atmospheres designed to re- 
produce the detailed chemistry, opacities, and emergent 
spectrum (but not the atmospheric motions), the effects 
of energy redistribution have been folded into a single 
parameter, referred to as a in this work. The a parame- 
ter is simply the ratio of the planet's cross-sectional area 
("TT-Rp) and the surface area of the planet from which the 
absorbed stellar luminosity is to be re-emitted. In this 
case, the stellar flux incident at the top of the model 
atmosphere becomes, 

Fin C ,x=a(^\ F+\. (5) 

For a more detailed description of the energy balance in 
an irradiated binary comp anion and t he de velopment of 
a similar a parameter, see iPaczvhskil ()1980l) . 

When the incident flux is scaled by a, the underly- 
ing assumption is that dynamical processes in the at- 
mosphere are efficient enough to uniformly distribute 
the incident luminosity over either the day hemisphere 
(a = 0.5) or the entire sphere (a = 0.25). In addi- 
tion, every point on the day side (or entire surface) is 
assumed to be identical and, thus, can be described by 
a single 1-D model with the same emergent and incident 
flux. Note that the a — 0.5 case corresponds to the aver- 
age (over /j,) of the incident flux defined in Eq. |2 Also, 
the no-redistribution and a = 0.5 cases receive and re- 
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Fig. 3. — Planet-star flux density ratios for day side models 
assuming no redistribution (top curve) and redistribution models 
with a = 0.5 (middle solid curve) and a = 0.25 (bottom solid 
curve). The model T-P profiles are shown in Fig. |5] IRAC band 
fluxes for each model (found by convolving with the IRAC response 
curves) are indicated with open squares and filled circles show the 
Spitzer data with lc error bars. The 4.5 /jm IRAC value for a 
10 X solar, a = 0.5, model is also shown (solid square). The lower 
dotted line corresponds to an isolated brown dwarf model with 
T cff = 1150 K. 



radiate the same amount of incident luminosity from the 
day side, but will predict very different T-P profiles and 
phase-dependent spectra (see below). 

3. RESULTS 

The greatest number of observational constraints exists 
for the planets TrES-1 and HD209458b. Models specifi- 
cally tailored for these two objects are presented below, 
adopting the parameters listed in Tabled an d are com- 
pared to the recent Spitzer data. 

3.1. TrES-1 

Figure [21 shows the sequence of solar metallicity T-P 
profiles across the day and night sides of TrES-1. On 
the day side, Ti, lt ~ 100K which is motivated by evo- 
lution calc ulations that reprod uce the observed radius 
of TrES-1 <|Baraffe et alJl2005|) . The night side atmo- 
sphere model with the same adiabat as the day side, has 
T c ff = 225K based on the entropy matching criteria men- 
tioned above. The sphere in Fig. [2] shows the concentric 
regions around the sub-stellar point represented by each 
atmosphere model. 

As expected, a steep temperature gradient along an 
isobar, VTp = (dT/dfi)p, is present from the sub-stellar 
point to the terminator. At P > 0.01 bar, VTp is very 
large across most of the day side and increases dramat- 
ically as the temperatures drop off near the terminator 
(designated by the night side T-P profile). Unlike the 
deeper layers, the top most layers of the atmosphere re- 
ceive a steady supply of incident stellar flux, even when 
a > 60° - however at shallow angles with respect to 
the surface normal. Consequently, VTp is smaller for 
P < 0.01 bar compared to higher pressure depths for 
most of the day side except near the terminator. Note 
also that the T-P profiles (for large a) become flatter 
at low P and steeper at high P and, eventually, become 
inverted very near the terminator. 
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Fig. 4. — Same as Fig. [2] but adopting the parameters for 
HD209458b in Table [I] 
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Fig. 5. — Planet-star flux density ratio for HD209458b assuming 
no redistribution (top solid curve) and redistribution models with 
a = 0.5 (top dotted curve) and a = 0.25 (bottom solid curve). 
The T-P profiles are shown in Fig. [I] Open squares indicate the 
predicted IRAC and MIPS flux densities. The observed 24 fim 
MIPS value is shown with la error bars (filled circle) along with 
the 24 (im MIPS response curve. The lower dotted line is the flux 
density ratio for a Brown Dwarf with T e ff = 1450K. 



The stellar heating also causes the atmospheric exten- 
sion 3 to increase from the terminator to the sub-stellar 
point by a factor of ~ 5. The change in extension coin- 
cides with a significant pressure gradient below the pho- 
tosphere along constant- height surfaces where the change 
in pressure from the day to night side can be factors as 
large as 50 to 100. Above the photosphere, the day- 
to-night pressure gradient is present, but fairly small. 
Despite the large increase in extension, the change in ra- 
dius is modest leading to no more than a 5% increase in 
the area of an isobaric surface on the day side compared 
to the night. Also of interest is the radiative-convective 
boundary which is not on an isobar and is significantly 
deeper at the sub-stellar point (P ~ 10 3 bar) compared 
to the terminator and night side (P ~ 10 bar). 

The condensation curves for Fe, Mg 2 Si04, MgSiC>3 and 
Na 2 S are also shown in Fig. [21 Cloud formation is typ- 
ically believed to occur near the intersection of the T-P 
profile and the condensation curve for a given species. 
While the models presented here are cloud-free, the con- 
densation curves suggest that photospheric clouds might 
be possible at a variety of atmospheric depths and com- 
positions but, in the no-redistribution case, would be 
confined to ~ 50% of the day side around the sub-stellar 
point. The no- redistribution model also suggests that 
clouds might form at much greater heights (i.e. lower 
P) around the sub-stellar point compared to the a = 0.5 
redistribution model. The presence of clouds would have 
an impact on the predicted T-P profile and, therefore, 
needs to be treated self-consistently. Cloud formation 
across the day and night sides will be explored in a fu- 
ture paper. 

The high day side temperatures lead to an atmospheric 
chemistry dominated by H2, He, H2O and CO. However, 
near the terminator and on the night side, most of the 
carbon is bound in CH4. In Fig. the planet-star flux 

3 The atmospheric extension is defined here as the difference 
between the radius at P = 10 — 6 bar and at P = 100 bar. 



density ratios 4 , 

l* UJ jv 1 j 

for the no-redistribution model are compared to the Sp- 
tizer observations along with standard 1-D models with 
a = 0.5 and a = 0.25. Since the stellar spectrum at IR 
wavelengths is fairly smooth, all of the features seen are 
due to absorption in the planet's day side photosphere. 
Strong absorption by H2O and CO are easily identified 
in the model. The two hottest models (no redistribu- 
tion and a — 0.5) are in good agreement with the obser- 
vations at 8 /im but significantly overestimate the flux 
density at 4.5 /im. The models with complete redistribu- 
tion (a = 0.25) agrees reasonably well at 4.5 /im and, at 
the 2a level, agrees with the 8 /im observations. Despite 
the broad wavelength span of the IRAC instrument, all 
of the IRAC bands probe a fairly narrow region of the 
planetary atmosphere between P = 10 _1 and 10~ 2 bar. 

Also shown in Fig. are the flux density ratios using 
a spectrum from a non-irradiated brown dwarf model at- 
mosphere with T c ff = 1150K which roughly corresponds 
to the equilibrium T e ff of the a = 0.25 model. While 
models of brown dwarfs appear to agree ver y well with re- 
cent S pitzer IRS spectra of L and T dwarfs ( Roellig et al. 
2004), the observations of TrES-1 are clearly inconsis- 
tent with a standard brown dwarf spectrum, especially 
around 8 /im. This disagreement confirms that irradi- 
ated EGPs have atmospheric structures very different 
from isolated brown dwarfs 5 ; a property that was not 
immediately obvious in the past. 

4 The ratio (Hp/H*) 2 , determined from the transit light curves, 
has error bars much smaller than the individual planet or stellar 
radii. Therefore, the flux density ratio comparison is not limited 
by t he typical uncertaint ies of an absolute flux determination. 

5 IFortnev~c t al. ( 200^) have also commented on the less than 
satisfactory agreement between models and the EGP Spitzer data 
in light of the very good agreement between models and Spitzer 
IRS brown dwarf spectra. 
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As demonstrated by our own solar system planets, 
it is possible that EGPs do not have the same metal 
abundances (or relative proportions) as their parent star. 
The main discrepancy between the no-redistribution and 
a = 0.5, solar abundance, models is the 4.5 /im flux 
which coincides with the strong fundamental CO absorp- 
tion band. The flux density ratio could be low at these 
wavelengths due to enhanced metal abundances. In or- 
der to achieve a CO absorption feature with flux density 
ratios as low as in the a = 0.25 solar abundance model 
(which reproduces the 4.5 /im IRAC observations), the 
metal abundance would have to be 10 times that of the 
parent star (see Fig. EJ). However, given the width of the 
IRAC band passes, the fluxes outside the CO absorption 
band remain high enough to keep the model's integrated 
flux density ratio above the 2a error bar. Note that in- 
creasing the C to O ratio does not improve the compari- 
son to observations. A larger C to O ratio increases CO 
absorption, but simultaneously lowers the water concen- 
tration thereby increasing the planet's fluxes at wave- 
lengths red-ward of the 4.5 /im that are also included in 
the IRAC band. 

3.2. HD209458b 

FigureEJshows the T-P profiles for HD209458b assum- 
ing solar abundances, and the parameters listed in Table 
^ The adopted inner adiabat was based on evolution cal- 
culations which suggest a substantial intrinsic luminos- 
ity for HD209458b' s mass and abnormally large radius 
(|Baraffe et al.ll2003D . The sub-stellar point T int = 230K 
and the entropy-matching non-irradiated model for the 
terminator and night side has T e g = 500K. 

The predicted T-P trend across the day side is simi- 
lar to that for TrES-1. The major differences are due 
to the greater parent star luminosity for HD209458b, 
which leads to a significantly hotter sub-stellar point. 
The intrinsic luminosity is also higher for HD209458b, 
which leads to warmer regions near the terminator. 
The decline in temperature above the nearly isother- 
mal photosphere (P < 0.01 bar) for the hottest por- 
tions of the d ayside is consistent with recent works by 
other groups llSudarskv et al.ll2003t iFortnev et al .1120051 
llro et al.l2005|) . This temperature decline points to some 
of the differences mentioned above between the present 
Rainout models and those based on the earlier Cond 
approximation. For HD209458b, with a = 0.5, the 
Cond assumption leads to a nearly isothermal profile (T 
~ 1700K) for most o f the atmosphere (see Fig. 2 of 
iChabrier et alJ l)2004|) ). While the complete removal 
of TiO and VO via the adopted Rainout process con- 
tributes to the cooler outer atmosphere and hotter pho- 
tosphere shown in Fig. 0] A more detailed comparison 
between Rainout and Cond models is given in Barman 
et al. (2005, in preparation). 

If clouds are sustainable on the planet's day side, then 
Fe, MgSiOs and Mg2Si04 clouds might form at very low 
pressures across most of the day side surrounding the 
sub-stellar point. Near the terminator, the temperatures 
are much cooler than the average temperature across the 
day side. Consequently, most of the limb is well be- 
low the condens ation temperatur e of Na, consistent with 
the findings bv llro et alJ lf2005). Condensation of Na 
near the limb may contribute to the lower than expected 
Na absorption detected with HST (Cha rbonneau et alJ 
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Fig. 6. — Planet-star flux density ratios for HD209458b assuming 
no redistribution at five different orbital phases from 0.44 (top 
curve) to 0.06 (bottom curve) and orbital inclination of 90°. Phase 
= coincides with the time of transit and only the night side is 
visible. At phase = 0.5 the planet is behind the star. 

l200l . 

The 24 /im MIPS observations probe the Rayleigh- 
Jeans tail of the planet's spectrum, in a region domi- 
nated by H2O line opacity. Figure [5] shows the day side 
planet-star flux density ratios for the no-redistribution 
model and 1-D models with a = 0.25 and a = 0.50. At 
the la level, the MIPS observations favor a strong re- 
distribution of the absorbed stellar flux (a = 0.25). The 
two hottest cases, no redistribution and a = 0.5, are only 
marginally in agreement with the 2a observational error 
bars. The MIPS observation probes the atmosphere at 
a pressure (P ~ 10~ 2 ) similar to those probed by the 
IRAC observations. Note also that for A > 10 /j,m, the 
planet's day side spectrum is nearly identical to that of 
a non-irradiated brown dwarf with T c ff = 1450K (which 
matches the emergent flux for the a = 0.25 case). The 
brown dwarf-like appearance of the spectrum at far-IR 
wavelengths is to be expected since, at these tempera- 
tures, 24 /im is well within the Rayleigh- Jeans tail. 

Another potentially useful limit for HD209458b has 
been s et at shorter wavelengths by ground based o bserva- 
tions (|Richardson et alJl2003HSeager et al.ll2005ft . These 
observations indicate the planet's spectrum may have a 
less prominent 2.2 /im peak (or an overall lower lumi- 
nosity) than predicted by many models. This limit also 
favors a < 0.50. 

3.3. Phase-dependent Flux densities 

Since the entire day hemisphere has been modeled by a 
collection of T-P profiles and intensity spectra, estimates 
for the phase-dependent fluxes can be constructed sim- 
ply by changing the star-planet-observer orientation and 
re-integrating the surface intensities (see Eqs. OHand0J|. 
Figure IH1 shows the predicted phase-dependent flux den- 
sity ratios for HD209458b with no energy redistribution. 
Note the significant drop in flux between 5 and 10 /im 
and the shift in the peak flux between 1 and 8 /im toward 
redder wavelengths. As the phase approaches zero, the 
coolest parts of the planet, which are dominated by CH4 
absorption, come into view. 

Flux densities for IRAC and 24 fim MIPS bands are 
shown in Fig. [7|for all phases with (a — 0.5 and 0.25) and 
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without energy redistribution. As the level of redistribu- 
tion increases, the thermal surface brightness becomes 
more uniform and the IR light-curves flatten out. For 
a = 0.25, the IR light-curves are constant with values 
equal to those shown if Fig. Note that optical and 
near-IR light-curves will not be flat when a = 0.25 due 
to reflected star light. 

A phase shift of TrES-l's hot side would not improve 
the simultaneous fitting of no-redistribution and a = 0.5 
models to the la observed 8 and 4.5 /im fluxes for 
TrES-1. For HD209458b, both a = 0.25 and 0.5 mod- 
els agree with the 2a error bars but favor a = 0.25 
at the la level. However, the large 2a error bars do 
leave room for a shift of up to 0.25 in phase - simi- 
lar to recent results from at mospheric circulation models 
ijCooper fc Showmanll2005|) . For both planets, the best 
fitting model appears to be one with a = 0.25. This 
may indicate that fast photospheric winds are present 
with speeds in excess of I km s -1 as predicted for 
atmos pheres with similar day-nig ht temperature differ- 
ences ijCooper fc Showm an 2005). However, using ap- 
proximations for the radiat i ve and advective ti me-scales 
i|Showman fc Guillotl 12002 ISeaeer et all I2005J) and as- 
suming 1 km s -1 wind speeds, T ra d/T a d v < 1 for P < lbar 
at the sub-stellar point for both TrES-1 and HD209458b. 
Small values for r ra d suggest that winds will not be capa- 
ble of entirely removing the large day-night temperature 
differences even at photospheric depths. 

4. DISCUSSION AND CONCLUSIONS 

Detailed radiative-convective equilibrium models have 
been presented above for the atmospheres of HD209458b 
and TrES-1. In the absence of energy redistribution, the 
models predict steep horizontal temperature gradients 
from the sub-stellar point to the night side that vary 
substantially with depth. Compared to earlier works, 
which assumed very efficient energy redistribution, the 
no-redistribution models predict much hotter tempera- 
tures across most of the day hemisphere and significantly 
cooler temperatures near the terminator. 

The existence of steep horizontal temperature gradi- 
ents in equilibrium models strengthens the case for strong 
zonal winds. However, fast winds (~ 2 km sec -1 ) ca- 
pable of altering the global atmospheric temperature 
profile are unlikely to entirely remove these tempera- 
ture gradients, especially at pressures lowe r than a bar 
l|Cooper fc Showman! 120051 llro et al J 12005(1 . bmce T ra d 
decreases rapidly with decreasing pressure (i.e. towards 
the top of the atmosphere), the horizontal and vertical 
temperature gradients may well have a mixture of the T- 
P characteristics shown above for the various redistribu- 
tion scenarios. For example, the temperature structures 
could be close to the static, no-redistribution, scenario at 
P < 0.01 bar where r ra d < 10 4 sec. Deeper into the at- 
mosphere (where r ra d >> r a d v ) advection is probably an 
important mechanism for energy transport which could 
lead to a T-P structure similar to that of an a = 0.5 
or a = 0.25 redistribution model for P > 1 bar. Conse- 
quently, the T-P profile could, in some cases, be fairly flat 
across much of the day side photosphere. Note that a flat 



photospheric T-P profile would produce a spectrum close 
to that of a blackbody - a possibility which, so far, is not 
exclud ed by Spitzer observations fsee lCharbonneau et alJ 
(2005) for a blackbody comparison to the TrES-1 Spitzer 
data). 

The 24 (Ltm Spitzer observations of HD209458b have 
been shown to favor an atmosphere undergoing efficient 
day-to-night energy redistribution. The 4.5 and 8.0 /im 
observations of TrES-1 also favor an atmosphere which 
is experiencing significant redistribution. However, an 
a = 0.25 redistribution model only agrees with both data 
points, simultaneously, at the 2a level. For both planets, 
it appears that the fully static, no-redistribution, case is 
ruled out by the Spitzer data. However, the data points 
are too few and have error bars that are too large to 
significantly constrain the model parameters. 

In anticipation of additional Spitzer observations, 
phase-dependent spectra have been calculated, along 
with phase-dependent planet-star flux density ratios. For 
HD209458b, placing limits on the fluxes at the quadra- 
ture and near night side phases (esp. at 24 /xm) would 
test the large night side luminosity (with T e g ~ 400 — 
500K) predicted by evolution calculations. The ratios of 
the planet's flux densities at phases 0.5 and 0.25 could 
also further constrain the degree to which energy is re- 
distributed to the night side. 

The current work is also applicable to planets in non- 
transiting orbits. There is a direct correspondence be- 
tween the phase-dependent fluxes shown above (for an 
edge-on orbit) and the inclination-dependent fluxes for 
a planet at superior conjunction in an arbitrarily in- 
clined orbit. Therefore, the predictions made above 
suggest that Spitzer could detect flux variations due 
to a close-in planet with orbital inclination as small as 
45° (i.e., phase 0.375 or 0.625 in Fig. Q. Since only 
two of the known transiting planets orbit stars bright 
enough to measure the planet-sta r flux density ratios , 
per forming similar obse r vation s as lDeming et al.l l)2005fl 
and lCharbonneau et alJ (|2005f) for non-transiting planets 
orbiting nearby stars would be very helpful. 
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